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Abstract 


In this research, we aim to analyze a mathematical model of Maize streak 
virus disease as a problem of fractional optimal control. For dynamical 
analysis, the boundedness and uniqueness of solutions have been investi- 
gated and proven. Also, the basic reproduction number is obtained, and 
local stability conditions are given for the equilibrium points of the model. 
Then, an optimal control strategy is proposed for the purpose of examining 
the best strategy to fight the maize streak disease. We solve the fractional 
optimal control problem by a forward-backward sweep iterative algorithm. 
In this algorithm, the state variable is obtained in a forward and co-state 
variable by a backward method where an explicit Runge-Kutta method is 
used to solve differential equations arising from fractional optimal control 
problems. Some comparative results are presented in order to verify the 
model and show the efficacy of the fractional optimal control treatments. 
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1 Introduction 


Maize is an important annual cereal crop of the world belonging to the fam- 
ily Poaceae. It is considered a staple food in many parts of the world. It 
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is the third leading crop in the world after rice and wheat [20]. Due to its 
highest yield potential among cereals, it is known globally as the queen of 
cereals. Maize streak disease (MSD) is the most serious viral crop disease in 
Sub-Saharan Africa. This disease is caused by the Maize streak virus (MSV), 
which was first described by the South African entomologist Claude Fuller in 
1901 [12]. MSV is mainly transmitted by as many as six leafhopper species 
in the Genus Cicadulina, but some other leafhopper species are also able to 
transmit the virus. In addition to maize, this virus can infect over 80 other 
species in the Family Poaceae. Severe MSD manifests as pronounced, contin- 
uous parallel chlorotic streaks on leaves, with severe stunting of the affected 
plant and, usually, a failure to produce complete cobs or seeds. Erratic epi- 
demics have been occurring every 3-10 years, and the main damage caused 
is to plants younger than six weeks old [24]. 

In recent years, mathematical modeling has become a valuable tool to 
study the mechanisms of plant disease spread, predict the future course of an 
outbreak, and appraise strategies to control. In most cases, differential equa- 
tions of the integer order have been used to construct such models; see, for 
example, [25, 8, 9, 14] and the references therein. The integer-order deriva- 
tives and integrals have local properties; that is, the next state is not influ- 
enced by the current and previous state. So, the integer-order mathematical 
models can not describe natural phenomena precisely. 

Fractional calculus is an extension of classical calculus that introduces 
derivatives and integrals of fractional order. Fractional derivatives have non- 
local properties, that is the next state depends on the current state and all 
previous states. This is the main excellence of fractional derivatives over 
classical derivatives. Due to this advantage, many applications of fractional 
calculus can be found in various fields of research, such as biology, economy, 
physics, control theory, and so on [15, 13, 21, 22, 26, 2]. In [23], a fractional 
model of tuberculosis disease has presented, and the values of parameters have 
been evaluated according to the actual clinical cases. In 2020, the dynamics 
of the fractional HIV infection model were studied by Evirgen Evirgen, Ucar, 
and Ozdemir [11]. Bozkurt et al., in their work [6], have analyzed a fractional 
model of COVID-19 by considering the fear effects of the media and social 
networks. In [5], the authors presented a fractional model for the simulation 
of the Cholera outbreak in Yemen. The authors in [4] proposed a fractional 
model to study the dynamics of the MSV in the maize plant population by 
considering the interaction of MSV pathogen with the past invasion. 

In light of this significant advantage, we were motivated to develop the 
model investigated in [3] into a new fractional model involving the Caputo 
derivatives. The Caputo derivative is of use for modeling phenomena that 
take account of interactions within the past and also problems with nonlocal 
properties. In this sense, one can think of the equation as having “memory.” 
After that, we discuss some properties of the fractional version of the model 
under consideration. Next, fractional optimal control (FOC) is applied as a 
generalization of the classical optimal control system [3]. The FOC model is 
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developed with three time-dependent control strategies proposed by Alem- 
neh, Kassa, and Godana [3]. 

The paper is organized as follows. In section 2, we give a brief review 
of the Caputo operator and discuss its basic characteristics. In Section 3, 
the fractional-order model formulation is presented, and the main properties 
of the fractional model are then given in section 4. Section 5 focuses on 
the dynamic analysis of the model. The FOC of the model and numerical 
simulations of the fractional model are presented in section 6. Section 7 also 
contains concluding remarks. 


2 Basic definitions and facts 


In this section, we give a brief review of the Caputo operator and discuss its 
basic characteristics [16]. 


Definition 1. For a function f : [0,t¢] ~ R,v € (n—1,n), and n EN, the 
left- and the right-sided Caputo fractional derivatives of order v of a function 
f are defined in the following forms: 


0D f(t) = iWon fe —u) "DF (u)du, t>0, (1) 
and ; 
— 1)" t 

£, Dz F(t) = = 1) 77 i (t—u)™ "DFM (udu, t<ty. (2) 


Here, I'(-) denotes the Gamma function. 


Definition 2. The integral operators related to (1) and (2), are specified by 


Cov F(t) = aC : (¢ — ny” (nan, (3) 
C5v F(t) = rw) zt "(9 — #24 (nan. (4) 


Additionally, if f € C™[a, 6], then 


6 to DFO) = FQ) —- ae (5) 
k=0 
nal (k) 
Pay, P04, F) = £0 — (ay a, 9 (6) 
k=0 : 
For any a1, a2 € Rand fi, fo € H'(0,ty), we have 
0 Df (ar fi(t) + a2 fo(t)) = a1 9 DF filt) + a2 F DY fro(t), (7) 
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Di, (arfi(t) + a2fo(t)) = a1 7 Dz, filt) + a2 7 Df, folt), (8) 
0 Jt (arfi(t) + aafa(t)) = a1 9 HA) + a2 6 Hf fold), (9) 
tt, (ar filt) + asfalt)) = an PI, filt) + 02 7H, folt). (10) 


Let f(t) be a constant function. Then 
0 Df (t) = 7 Di, f(t) =0. (11) 


The Caputo derivatives satisfy the Lipschitz condition. 


3 New fractional model of MSV disease in maize plant 


In this section, we develop a deterministic eco-epidemiological fractional 
model for the dynamics of MSV disease in maize plants. The original version 
of this model is a system of ordinary differential equations that have been be- 
fore presented in [3]. The effect of previous states in the current states of the 
disease spread has not been considered in this model. One way to overcome 
this drawback is to replace the integer-order derivatives in the model with 
noninteger-order derivatives [19]. Hence, we replace the ordinary derivative 
with the following Caputo fractional derivative operator 


d lo 

di piv Di, (12) 
where the auxiliary parameter 9 > 0 represents the fractional time compo- 
nents in the system. Thus, the new model is described by the system 


0’ CD" S(t) = rS(1 ae) a 
mgar itt) = FS — jal, a 
oS DY H(t) = q-— BF — wo, 
oF DtY (t) = on usY, 
S(0) = So, (0) = Io, H(0) = Ho, Y(0) = Yo, (14) 


where 0 <a < 1, Ni (t) = S(t)+J(t), No(t) = A(t)+ Y(t), and (S,/,H,Y) € 
R4.. In this model, S(t) denotes the density of the susceptible maize, and 
I(t) denotes the density of the infected maize. The susceptible and infected 
leafhopper vector densities are denoted by H(t) and Y(t), respectively. All 
parameters in the model are nonnegative. Description of the parameters are 
found in Table 1. 
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Table 1: Explanation of MSV model parameters 


Parameter Explanation 
By Predation and infection rate of infected leafhopper on susceptible maize plant 
Bo Predation and infection rate of susceptible leafhopper on infected maize plant 
A The half-saturation rate of susceptible maize with infected plant 
C The half-saturation rate of susceptible leafhopper with infected maize plant 
Kk Carrying capacity 
q Recruitment rate of susceptible leafhopper 
b Infected leafhopper conversion rate 
r Maize population intrinsic growth rate 
LM Death rate of infected maize 
[2 Death rate of susceptible leafhopper 
Ls Infected leafhopper death rate 


As can be observed, model (13) involves a system of nonlinear fractional 
differential equations. The exact solution of this model may not be available 
in general. However, a mathematical analysis of the existence and uniqueness 
of the solution ensures that a unique solution exists under some conditions. 


4 Properties of the model 


In the following, the model’s main properties are provided. Our model can 
be formulated as 


6 Di V(t) = O(t, V(t), 
V(0) =Vo, (15) 


where V(t) = (S(t), I(t), H(t), Y(t)). 
Lemma 1. [17] Let w(t) be a continuous function on [tp, 00) and satisfying 


6 Diw(t) < —Aw(E) + p, 


w(to) = wo, 06) 


where 0 < v < 1,(A,) € R?, A #0, and to > 0 is the initial time. Then 


Lb 


w(t) $ (wo ~ 5) Bu[-AE— to)"] + 5 


Xd ’ 
where E,, represents Mittag—Leffler function. 


Lemma 2. [7] Let 0 <y < land A< 0. Then £,,,,(At”) tends monotonically 
to zero as t + oo. 


Lemma 3. [10] Let ® : [tg,0o) x R” > R” be a continuous function and 
Lipschitz-continuous respecting to the second variable. In addition to, let 
v € (0,1) and Vo € R”. Then, the problem 
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COYV(t) = B(t,V(t)), £> to, 
V(to) = Vo, (17) 


has a unique solution in C'([0, 00); R”). 


Theorem 1. All solutions of system (13) that initiate in R4 are bounded 
within the region defined by 


re 
OQ = {(S,I,H,Y) € R4|S(t) + I(t) + H(t) + ¢$Y() < : +e, for all e > O}. 


Proof. Define a time-dependent function w(t) = S(t) + I(t) + H(t) + ¢Y(t). 
So, for any positive number p, we have 


0 Dy w(t) + pw(t) 


=rs(1 ott) jal — po PY + pS+ pl +pH+FY +q 
<rS(l *) t(p—pi)l+(p—p2)H+(p Hs),¥ +g 

= (r +p) — 28? + (p—p)E + (p— wa) + (9 13) ,¥ +4 
S a tp)” + (p— pa) + (p — pa) + (p Hs);¥ + g. 


Taking p < min(pu1, Ja, [43), So 
Cav w(t) + pw(t) < L, 


where L = f(r +p)? + q. Now, we apply Lemma 1 and obtain 


w(t) < (w(0) — 2)E,[-pt”] + 2 


L 
Thus, w(t) > — as t > oo and 0 < w(t) < —. Hence all solutions of system 


(13) that starts from R4 are confined in the region Q = {(S,I,H,Y) € 


L 
R4 |w(t) < r +e, for all ¢ > O}. 
Now, we study the existence and uniqueness of system (13) in the region 
A x [0,T], where 
A = {(S,1, H,Y)eR® : max(|$], |Z], |H|,|Y¥1) < M}, 
T < oo and M is sufficiently large. 


Theorem 2. For any nonnegative initial conditions, system (13) has a unique 
solution. 
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Proof. Let X = ($,J,H,Y). Consider a mapping 


Q(X) = (Qi1(X), Q2(X), Q3(X), Qa(X)), 


where 
Qu(x) = ng — EA) BE, 
Qo(X) = FEE — yal, 
IH 
et ea aft 
For any X,X € A, we have 
ax) - OI 
= |Qu(X) — Q1(K)| + 1Q2X) — Qa] + 1QaX) — Qa] + 1QalX) - Qa] 
=|rsq- 544) ae rB(1 oat) 4 ASE 
=|r(s -5)- 2s? 5%) - Lsr-3D) - Gg - 
you ey - 22) - mt - D4 - a PG - A) - at - 
opa( A — F ) _ sy -¥) 
< (r+ M53) + 28g - +m tole - oe 
+ pall — I| + p2|H — A\+ psl¥ — | 
Bt SM 28M g cae ey) 
(M4 6) tyne) + (ee 6) 4 a] 
< HIIX — Xl, 


where 


3rM_—-.28,M  26;(A+M)M 


H = max{r 4 K A ’ A2 + ps3, 
M(1+6 M(C+M)(1+ 6 
aMA+8) , ,, BMICHMO+ | 
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Thus, Q(X) satisfies the Lipschitz condition with respect to X. Hence 
there exists a unique solution of the system (13) with conditions (14) on 
A x (0,7). 


5 Dynamical behaviors 


One of the key concepts in epidemiology is the basic reproduction number 
(BRN). The aim of this section is to obtain the BRN for model (13) and study 
the local stability behavior of the model at its disease-free equilibriums. 


5.1 Basic reproduction number 


Consider the following fractional differential system: 


o” 1§ D7 X(t) = F(X,Y), 
o” “FDL (t) = G(X, Y), (18) 
G(X,0) =0, 


with nonnegative initial conditions X(0) = Xo € R? and Y(0) = Yo € R?, 
where the components of vector X = (S, H) represent the number of suscep- 
tible maize and leafhopper, and the components of vector Y = (1, Y) indicate 
the number of infected maize and leafhopper. Furthermore, we presume that 
the function G is of class C!, F is continuous, and the system (18) with the 
initial conditions X (0) = Xo and Y(0) = Yo admits a unique solution. Also, 
suppose that E = (X*,0) € R* denotes the disease-free equilibrium point of 
the system (18). Let A = oG (X*, 0) = M — D, where M, D are two square 
matrices that D > 0 is a diagonal matrix and M > 0. Then the BRN Ro is 
obtained as the spectral radius of MD~!. 

For system (13), we have 


S BS 
A= | spt 3] =| po -|¢ 2 (19) 
(C+? —H3 (C+? 0 0 3 


So, 


0 CAS 81 82HCS 
R = p3(A+S) = 1P2 . 5) 
a (ie 0 pip3(A + S)(C + I)? (20) 
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5.2 Local stability analysis 


Theorem 3. The disease free-equilibrium point Eo = (0,0, ar) of sys- 
tem (13) is always unstable while the disease-free equilibrium point Ey, = 
(K,0, mE 0), is locally asymptotically stable if Ro < 1. 

Proof. Conforming to Mittag-Leffler function [18], the disease free equilib- 
rium F of system (13) is locally asymptotically stable if all eigenvalues 
Ai,t = 1,2,3,4 of Je satisfy jarg(Ai)| > 41 = 1,2,3,4. The Jacobian 
matrix associated to Eo is given by 


r O 0 O 
=< @ o 

JE = 0 —B24q =P) 0 
0 


The eigenvalues of the matrix Jz, are Ay =r >0, Ag =p <0 A= 
—p2 <0 Ag = —p3 < 0. We observed that |arg(\1)| = 0 < 4. So, the 
equilibrium point Ep is unstable. 

The Jacobian matrix associated to EF is the following one: 


—-BK 

= = 0 ALK 

0 -— go 6A. 

JE, = o ae 
1 0 c = —pl2 0 

0 GS 0 -ps 


The following characteristic equation is obtained from Jz,: 


60) = (r+ Nee +A) (28+ Ge ts) + Gets — eet) 


We observe that two roots of the characteristic equation $(A) are 


M=—r<0, Ae=—pe <0. 


It is obvious that |arg(A1)| > 4 and |arg(A2)| > 4. The remaining eigen- 
values are given by 


Kbqhi Be 


»? ee pe 
+ (fa + M3 )A + (Hi pb3 (A+ Kyi 


(21) 


By the Routh—Hurwitz criteria, all the roots of the polynomial (21) are neg- 
ative or have negative real part if and only if 


_ _ Kbqh1 82 
LH1ip3 (A+ K)jisC ’ 
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or 


Kbq81B2 
Ro = <i. 22 
: la + K) ui pop3C (22) 


Hence, FE; is locally asymptotically stable if Ro < 1. 


6 Optimal control problem 


In this section, to attain the minimized number of infected maize and infected 
leafhoppers, we reconsider the model (13) and formulate an optimal control 
problem with three control variables u1(t), u2(t), and ug(t). Let 


U = {(u1, v2, us)|u1,u2, and uz are Lebesgue measurable on [0, ty], 
0 < wu, U2, u3 <1, for all t € [0,ty]}, 


be the admissible control set. With the existence of control u1, it is expected 
to diminish the number of infected maize as compared to those without con- 
trol cases. The control variable w2 is used to control the number of infected 
leafhoppers. Furthermore, u3 is chemical control that is used as an inter- 
vention strategy to optimize the objective functional *. After incorporating 
the control variables u(t), u2(t), and u3(t) in the model (13), the optimal 
control model is as follows: 


eee Ga.7 B, SY 
oS DYS(!) = rS(1 ae . Gta) ot 
oan = ars 3 (u1 + U2)L, 
o’ OD’ H(t) =q-(1 2) ear (ug + 2) H, (23) 
y-10@VvV (t) = (1 — ug) oat 7 7 (us + us)¥, 


(0), 1(0), H(0), ¥(0) > 0. 


Consider the following objective functional: 
i 1 2 2 2 
F= | djl +d2Y + 5 (wii + wous + w3u3) |} dt, (24) 
0 


where dj,d_ are the weights on the state variables and w,,w2, and w3 are 
relative weights of the treatment related to the control functions u,, ug, and 
U3. 


Our aim is to minimize the cost value F by the state and control variables 
I*,Y*,uj,u3, and uz satisfying the constraints (23). For this purpose, we 
use a kind of Pontryagin maximum principle in the fractional order state [1]. 
We define the Hamiltonian function as below: 
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1 
§(S,I,H,Y)=d\I+d2Y¥ + 5 (wits + woud + w3u3) 


+m(rs(.- 25°) (1 - mn) SS) 
+ eo((1 =) Fg = (i + ua) 
+ slg (1 ua) 2 = (us + 12) 
+ wa((1 ~ 2) BEE — (us + 15)¥)) 


where w;,i = 1,2,3,4 are the co-state variables or adjoint variables. The 
optimality conditions are obtained from 


09 OH OH 


— — = 0. 
Ou, Oug Ou3 
Hence, we have 
_ Bo(we = m1)IH 
wi(A + S) : 
wot Bo(bwa = w3)lH 
U2 = ; 
: we we(C + I) 
Y — w3 
lige (25) 
ws 
where the adjoint variables satisfy 
7 ay K—2S-1 BiYA 
Di, 71 (t) = aS ={r( K )-a “tay sym 
Bi(1 — u1)SY 
+t aay 
OS) Ss CH 
PE, mall = Gp me — Fem — in Fada — OG 
BobCH 
+ (1 — uz) (c+ pe 
(e) I bb — I 
PDE, walt) = OP ={(ua — 1) PE — na — us }ems + “Oy, 
0. —1)s 
P Di, walt) = =d2 4 att (wi — w2) — (ua + 3) wa, 
wi(tr) = wa(ty) =w3(tp) = wa(ty) = 0. (26) 


Then, we have the following boundary value problem for optimal treatment: 
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v-1C gy aed B, SY 
o 18 D¢S(t) =rS(1— =) -A-m) as 
By SY 


Ww 


oF D¢T(t) = (1 eee + us), 
of DYH(t)=q-(1 uy) B (ug + 2) HA, 


C+I 
b6,1H 
aay) = (1 — ug) m — (ug + w3)Y, 
, K-25 — AYA 
POY, a(t) = {r( K pei “AE 52 jo 
B01 — u,)SY 
Dt, walt) =d, — 2a — (m1 +u2)m2-(1 us) er W3 
1m Bie a 
u 
Dt, w(t) = {(u2 a fa — Ug}m3 + pg 
, By (ur = 1)S 


YDi,wa(t) = de 4 (m1 — we) — (us + 3) 4, 


A+S 
mr (ts) = wa(ty) = walts) = wal(ty) =0, 
S(0) = So, 1(0) = Io, H(0) = Ho, ¥(0) = Yo, 


where u(t), u2(t), and u3(t) are given by (25). In turn, the optimality condi- 
tions of Pontryagin’s Minimum Principle establish that the optimal controls 
uj(t),us(t), and u3(t) are defined by 


Bo(we2 = o1)1H 
~ wi(A+ sy fh 


uz; = min{max{0, 


> mol Bo(bwa = w3)lH 
= 0 1 
Ug min{max{ ? we = wo (C A I) Ts fy 
Y — 73H 
uz = min{max{0, ck EL 1}. 


W3 


Simulation and discussion 


In this part, the effects of fractional operators on the behavior of controlled 
system for the dynamics of MSV disease are investigated. We develop the 
fractional version of fourth-order Runge-Kutta (RK4) algorithm for the cou- 
pled system (27) and apply the iterative process as follows: 

We use S(0) = 1000,/(0) = 20,H(0) = 100, and Y(0) = 0 as initial 
values. In addition, the parameter values can be seen in [3]. 

The dynamical behaviors of all variables in the new fractional model with- 
out applying any control for different values of the fractional orders and the 
classic integer-order are plotted in Figure 1. As seen in this figure, infectious 
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Algorithm 1 


Step 1 Set the initial values for the control functions wuj(t), u2(t), and u3(t). 

Step 2 Use the current values of control functions and apply the forward 
fractional RK4 method for the control system and obtain the original 
variables. 

Step 3 Apply the backward fractional RK4 method to compute the adjoint 
variables using the current values of the original variables and control 
functions. 

Step 4 Update the value of control functions. 

Step 5 If the updated values of the original variables, adjoint variables, and 
control functions are not close enough to their previous values, then go 
to Step 2. 


maize and leafhopper densities increase with the fractional orders decrease 
and tend uniformly to the integer-order trajectory. Furthermore, when the 
fractional orders decrease, the densities of susceptible maize and leafhopper 
are reduced and go to the v = 1 state. 

To indicate the efficiency of the new optimal control model, the same 
impact rate has been considered for all three controls, and the numerical 
results of the new model are compared with the classical integer model, in 
Figure 2. As can be seen in this figure, the participation of controls leads to 
a further reduction of infected maize in the new model than in the classical 
model. Therefore, the effect of controls on the fractional system is more 
successful than applying controls on the integer system, and the difference 
between them is significant. Of course, it should be noted that if no control 
is applied, the fractional model still leads to a significantly lower infection 
density than the integer model (Figure 2). 

In the following, we numerically examine the effect of several optimal 
control scenarios, where each scenario includes more than one interventionist: 


Scenario 1 Applying quarantine (wz) and chemical control (w3) along with 
elimination of prevention (u1). 


Scenario 2 Applying prevention (u;) and quarantine (wz) along with elim- 
ination of chemical control (u3). 


Scenario 3 Applying prevention (u,) and chemical control (u3) along with 
elimination of quarantine (wu). 


Scenario 4 Applying all three controls u,, u2, and us. 


In the first scenario, the prevention (ui(t)) effect is removed, and two con- 
trol functions u2(t) and u3(t) are used. Figure 3 shows that in this control 
scenario, the number of infected maize decreases, while if no control is ap- 
plied, the number of them increases over time. In addition, for fractional 
derivatives with lower orders, the rate of reduction of infectious cases is more 
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Infectious Leafhopper 
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2 f 2k 6 1 en 18 
Time(months) ‘Time(months) 


Figure 1: Numerical solutions for classical and fractional order models without controls 


significant. Hence, this scenario is effective in the decline of infection in the 
maize community, especially in the fractional model. 


The results of scenario 2 are presented in Figure 4. In this case, the 
quarantine control (u2(t)) is maintained as in the previous scenario, but the 
chemical control (u3(t)) is replaced by the prevention (u(t)). This scenario 
prevents the spread of infected maize and reduces their number. Therefore, 
this strategy is also successful in eliminating the disease in the maize com- 
munity. 

Figure 5 indicates the results of using scenario 3. In this scenario, the 
quarantine control is removed, unlike the previous two scenarios, and the 
rest of the controls are applied. Based on this figure, the number of infected 
is reduced compared to the case where there is no control, and this scenario 
also limits the growth of infectious cases. 


Finally, all control interventions are considered together. As you can see 
in Figure 6, the use of these controls is a successful plan and causes the 
infected maize to be destroyed as passing the time. 
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Based on Figures 3-6, the quarantine control is more effective than other 
controls. Without quarantine control, despite preventive and chemical con- 
trols, the number of infected maize will spread, but this increasing process is 
much slower than the case where there is no control, and also, the difference 
between their increasing manner is very significant. 


Infectious Maize 


‘Time(months) 


Figure 2: Numerical solutions of I(t), with uncontrolled and controlled conditions for 
classical and fractional order models 


‘Time(months) 


Figure 3: Numerical solutions of /(t) in classic and fractional model, with quarantine 
and chemical controls (ug 4 0,u3 #0, ui = 0) 


‘Time(months) 


Figure 4: Numerical solutions of I(t) in classic and fractional model, with prevention 
and quarantine controls (ui 4 0,u2 4 0, ug = 0) 
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‘Time(months) 


Figure 5: Numerical solutions of I(t) in classic and fractional model, with prevention 
and chemical controls(ui 4 0,u3 4 0,u2 = 0) 


1 
‘Time(months) 


Figure 6: Numerical solutions of [(¢) in classic and fractional model, with prevention, 
quarantine and chemical controls(u; 4 0, ug 4 0,u3 4 0) 


7 Conclusion 


In the present study, we developed a new mathematical model involving 
the Caputo fractional derivative for MSV disease in maize plants. First, we 
proved that the solution of this model system exists uniquely and that all so- 
lutions remain positive and bounded whenever they start with positive initial 
values, thus justifying the well-posedness of a biological model. We also de- 
termined the BRN for the model. Then, we studied the local stability of the 
disease-free equilibrium points of the model. The study demonstrated that 
one of the equilibrium points is always unstable, and the other equilibrium 
point is locally asymptotically stable if the model’s BRN is less than unity. 
Next, an optimization problem is formulated. Our main focus in this work 
is to investigate the influence of fractional-order derivatives on the optimal 
control problem. The optimality system was solved numerically by use of a 
forward and backward RK4 scheme. The effectiveness of preventive, quar- 
antine, and chemical controls on the fractional model is investigated in the 
figures. Different scenarios for the participation of these controls were eval- 
uated for various fractional-order values. We observed that in all scenarios, 
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the efficiency of the controls increases by moving away from the integer-order 
and reducing the fractional orders. Moreover, it was observed that the quar- 
antine control is more effective than other controls. Without quarantine 
control, despite preventive and chemical controls, the infectious density of 
maize progresses with an increasing trend. Of course, it should be noted that 
the increasing trend is much slower than the case where there is no control, 
and also, the difference between their increasing manner is very significant. 
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